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Abstract 


A rigid object of general shape is fixed inside a wind tunnel. The drag force 
exerted on it by the wind is determined by a new method based on simple 
basic Physics concepts, provided one has a solver, any solver, for the corres¬ 
ponding dynamic Navier-Stokes equation which determines the wind velocity 
field around the object. The method is completely general, but here we apply 
it to the traditional problem of a long cylinder perpendicular to the wind. 
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1 Introduction 


Consider a rigid object facing the counter flux of a viscous fluid (air or water, 
for instance). The fluid velocity held around the object (which would be uniform in 
its absence) is distorted by its presence. This velocity held v(r) can be obtained as 
a function of time t by solving the so-called Navier-Stokes dynamic equations (1, 2). 
These equations can be written in different forms, here we adopt a simple one 

(1) 

see for instance (3), where 

fl = Vxi) (2) 

is the vorticity held. Also, provided all speeds are small compared with the sound 
propagation speed in the same fluid, it can be considered incompressible, he., 


V •v = 0. (3) 

The above equation (1) is already expressed in adimensional units, conside¬ 
ring the wind speed V far from the object and some characteristic linear dimension 
of the rigid object such as the diameter D of a cylinder, for instance. Together 
the characteristic fluid properties p. its uniform density, and /i, its viscosity, these 
parameters can be condensed into a single adimensional Reynolds number (4) 


Re 


pVD 

p 


( 4 ) 


This way, one can consider V in equation (1) as the speed unit, and D as the length 
unit. The corresponding time unit is D/V. 

These 4 above equations determine completely the velocity or vorticity holds 
around the object, as functions of time, given the proper boundary and initial con¬ 
ditions. Normally, the boundary condition is v — 1 far from the object along some 
hxed direction, and v = 0 at the object surface. The initial condition may be (as 
here) the wind tunnel initially switched oh and suddenly switched on with a hxed 
value V (or Re, impulsive switching). In this case, one has hrst to solve the time 
independent Stokes equation V 2 f2 = 0, obtained by taking Re —X 0 in equation (1), 
in order to obtain the initial helds v(r) and O(r) at t — 0. By neglecting this care 
with the initial helds, one obtains spurious transient helds just after t = 0, while 
the steady-state final regime is not reached (5). 

What is not explicit in the above equations is how to determine the drag 
force the wind exerts on the hxed object. The traditional way to do this is by hrst 
determining the pressure distribution and the shear stresses on the object surface 
and then integrating them. However, these distributions are often extremely difficult 
to obtain, either experimentally or theoretically (6). Alternatively, one can obtain 
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the drag force just through the knowledge of the fluid velocity held around the 
obstacle. Also in this case, an integration over the surface of the object that involves 
the gradient of this held on the surface is necessary (7, 8). If the method adopted to 
solve the above equations is a numerical one on a discrete grid, the determination 
of the quoted gradient requires a very fine grid near the surface, which takes a large 
computational effort. That is why normally researchers adopt non-uniform grids, 
very fine only near the object surface, when the aim is to determine the drag force. 
This approach unnecessarily complicates the mathematical discretization procedures 
needed to translate the continuous derivatives into finite differences, besides the 
further computer effort fine grids require. 

The current text introduces an alternative method to determine the drag 
force replacing the surface integration by a volume integration, which solves both 
problems commented at the end of last paragraph. The method is completely gene¬ 
ral and based only on simple basic Physics concepts. It is described below, and its 
validity mathematically demonstrated at the Appendix (although this mathemati¬ 
cal demonstration is completely unnecessary, the conceptual arguments below are 
enough). 

Consider one has a solver for the Navier Stokes equations above, providing 
the held v in all space at time t + 5t from the knowledge of the previous held at 
time t, where 5t is some small time interval. Any solver can be used, the accuracy 
of the method now described depends only on the accuracy provided by this solver. 
Imagine one replaces the rigid object by an extra portion of fluid at time t. This 
extra portion of fluid is also at rest in the same way as the removed object. But 
it would move from time t on. A small fluid velocity would appear at each point 
inside the volume formerly occupied by the object, at t + 5t. By integrating these 
velocities inside the quoted volume and multiplying the result by the density of the 
fluid, one obtains the total linear momentum which would be transferred from the 
fluid to the object, a vector. This would-be transferred momentum is compensated 
by the mechanical device keeping the object at rest inside the wind tunnel. Finally, 
dividing this total linear momentum by 5t, one obtains the force exerted by the 
original fluid on the object. 

In short: I) One has the velocity held at time t around the object, with all 
other velocities null inside the volume formerly occupied by the object, now occupied 
by the static fluid replacing it; II) One uses the Navier Stokes equations solver in 
order to obtain the same held at t + 5t] III) Non-null, small velocities appear inside 
the quoted volume; IV) One integrates these penetrating small velocities over the 
quoted volume, multiply the result by the fluid density and divide it by 5t. 

It is worth mention why the Newtonian momentum transfer procedure is 
enough in our case, no convective contributions are needed. There are two different 
descriptions to treat continuous fluid movements. The Lagrangian description takes 
the control volume of fluid moving with the fluid, the portion of fluid inside this 
volume is always the same, nothing entering or leaving this control volume. In this 
description, one can directly apply Newton’s law. The Eulerian description takes a 
hxed and static control volume, the fluid passing through it. The relation between 
these two descriptions is made through the so-called convective derivatives. The 
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material derivative (D/Dt) applicable to the Eulerian description is obtained by 
adding the quoted convective derivative (v x d/dx + v y d/dy + v z d/dz) to the regular 
time derivative ( d/dt ) applicable to the Lagrangian description (9). This extra term 
is just the origin for convective corrections over merely Newtonian inertial forces. 
However, in our particular case, by replacing the object by static fluid, our control 
volume is that of the object itself. Because the fluid is static at this precise time 
t, the Lagrangian and Eulerian volume controls are exactly the same, there are no 
translations between them {d/dx = d/dy = d/dz = 0), therefore the convective 
derivative of any quantity vanishes. That is why convective terms don’t enter in 
our momentum transfer formulation. This also applies to the case of a rotating long 
cylinder perpendicular to the wind treated later, the convective derivative is also 
null in this case. 

This text is organized as follows. In next section we show the results for the 
drag force on a long static cylinder, perpendicular to the wind, inside a wind tunnel 
initially switched off and suddenly switched on in t — 0 with a Reynolds number 
Re = 1, 000. After some transient time, the steady state drag force reaches a value to 
be compared with the experimentally known value of the drag coefficient Cd ~ 1.0 
(10, 11, 12, 13, 14, 15). The agreement with this experimental value is excellent 
(see figure 3 discussed later). Before that, however, starting at t — 0 the drag force 
increases and reaches larger values, decreasing a little bit up to stabilizing at the 
steady state. This behavior is the same observed in experiments as (16, 17, 18, 19, 
20, 21, 22), in particular that of a small and light falling ball which surprisingly 
reaches a larger speed than its final limit speed during the transient time (23) . The 
explanation of this curious phenomenon resides on the gradual formation of the von 
Karman street of vortices behind the ball: While it is not yet completely formed, the 
drag force is smaller than its own steady state value, allowing larger fall speeds; After 
the street of vortices is completely formed, the drag force surpasses the ball weight 
and the ball breaks; At the end, the drag force decreases a little bit becoming equal 
to the ball weight, and hereafter the limit speed remains constant. In the following 
section, a rotating cylinder is treated. Besides the drag, the lift coefficient due to the 
Magnus effect is also determined for various angular velocities of the cylinder. In all 
cases, the formation of vortices behind the cylinder is dynamically shown. Finally, 
the last section resumes the whole work. 
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2 Drag force on a long static cylinder 

In this section, one describes the flow over a long static cylinder immersed 
in a wind tunnel initially switched off and suddenly switched on in t — 0 with a 
Reynolds number Re = ROOO 1 and, then, one applies the method described in the 
previous section to determine the behavior of the drag force on this cylinder. 

Starting from the Stokes configuration, the wind is switched on at t — 0 
with a fixed Reynolds number Re = 1, 000. After a transient time, two symmetric 
vortices are formed behind the cylinder, a close up of one of them is shown in figure 
1. Some time later, they move a little bit downstream and start to become stretched 
along the X direction, as shown in figure 2. 



Figura 1 - Starting from the Stokes configuration, the wind is switched on at t — 0 
with a fixed Reynolds number Re = 1, 000. The figure show the result 
at t — 1, when a fluid element far from the cylinder has already traveled 
one diameter. Just a close up of the region behind the cylinder, above 
the X axis, is shown. The other symmetric vortex below is not shown. 
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Figura 2 — Later than figure 1, t — 2. Vortices are still symmetric. 


After that, one vortex bifurcates and becomes a pair of vortices running in 
the same sense and repelling each other along the X direction. Later on, so does the 

1 One uses the method of successive relaxations to solve the Navier-Stokes equations (24). 
However, any other method may be used. 
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other vortex. Only then the aforementioned von Karman street begins to be formed: 
one vortex running clockwise slowly goes away downstream, followed by another 
running counterclockwise, and so on. 

A completely description about the flow on a static cylinder can be found in 
reference (24). There, the velocity held is considered two dimensional on the plane 
perpendicular to the cylinder, a traditionally adopted good approximation for these 
not-so-high Reynolds numbers. The external boundary condition is a lOD x 5 D 
rectangle, measured in terms of the cylinder diameter D. Outside this rectangle the 
velocity held V is considered uniform along the direction of its two parallel largest 
edges. The cylinder center is located inside this rectangle at a distance 3 D from the 
incoming wind edge and 2.5 D from both largest edges. The discretization numerical 
grid covering this rectangle has 400 x 200 pixels, therefore each pixel is a small 
x A/} square. The discretized time interval adopted to solve the Navier-Stokes 
equations is 4 ^ 7 , the time wind spends to travel one tenth of a pixel. 

So, one may now apply the previously described method to determine the 
behavior of the drag force on the static cylinder. Figure 3 shows the drag force as a 
function of time . 2 



Figura 3 - Drag force on a static cylinder in the wind tunnel initially switched 
oh. It is turned on at t — 0 with a constant Reynolds number Re = 
1,000, which defines the speed of the wind. The force in the direction 
perpendicular to the wind how is negligible. The dimensionless drag force 
is plotted in terms of the so-called drag coefficient. During one time 
unit the wind travels one cylinder diameter. The force increases, reaches 
a maximum value, decreases and finally stabilizes, corresponding to a 
dynamic situation where successive vortices appear continuously, slowly 
going away, one running clockwise, the next running counterclockwise 
and so on. The long von Karman street is then formed. Inset shows the 
same in a longer time scale. 


The time St adopted is the same as the discrete time interval used to solve the Navier-Stokes 
equation, i.e. the time the wind takes to transverse 0.1 grid pixel. 
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Analysis of the graph shows that the drag force on the static cylinder in¬ 
creases, reaches a maximum value, decreases, and stabilizes. As expected, the same 
behavior of the falling ball (23) is found for this large Reynolds number, for which 
the long von Karrnan street of successive vortices is present. On the other hand, 
for small values of Re, about a few dozens, for which instead of the long street pe¬ 
riodically fed with new vortices only two vortices appear and stabilize behind the 
cylinder, the force does not present a maximum value before stabilization. Also, in 
this case the drag force is much smaller. 

The existence of a transient period during which the drag force reaches a 
value larger than its steady state is confirmed by various experimental studies (16, 17, 
18, 19, 20, 21, 22). Moreover, our results are in accordance with known experimental 
data. A comparison between the results obtained by our method and experimental 
data is made through the drag coefficient (6), 

^ , _ Rdrng 

D = \pV 2 DV 

where Fd ra g is the parallel-to-the-wind component of the drag force, p is the density 
of the fluid, V is the wind’s speed, D is the diameter of the cylinder and l is the 
length of the cylinder. Drag coefficients are measured experimentally since the early 
twentieth century. The drag coefficient measured at steady state for a static cylinder 
was measured for the first time (within the knowledge of the authors) in 1921 (10). 
After that, several other similar experiments were carried out, and today it is a 
known fact that for Re = 1,000, Cd ~ 1.0 (11, 12, 13, 14, 15), whereas our method 
provides Cn = 1.01. 


(5) 


3 Drag force on a long rotating cylinder 

In a different version of the problem, the cylinder rotates with constant angu¬ 
lar speed to. Now, the fluid velocities on the surface of the cylinder are not null, but 
are equal to the speed at the surface of the cylinder, i.e., the fluid rotates with the 
cylinder. The angular speed to is measured in units such that to = 1.0 corresponds 
to the magnitude of the velocities on the surface of the cylinder equal to the speed 
of the wind. 

Let’s start with the Stokes laminar limit (Re —>■ 0) for which the Navier- 
Stokes equation reduces to the Laplace equation, V 2 D = 0. The resulting streamlines 
are shown in figure 4. 

When the original Navier-Stokes equation is solved for a fixed and large 
Reynolds number, say Re = 1,000, the solution shows that, after a transient time, 
as in the case of the static cylinder, the system cyclically stabilizes: a vortex rotating 
in the counterclockwise direction appears in the region behind the cylinder and below 
the horizontal axis formed by its diameter, see figure 5(a); then, this vortex grows 
and slowly moves away downstream; after some time, a second vortex appears in 
the region above the horizontal axis, rotating in the clockwise direction, see figure 
5(b); this vortex also grows and moves away; later, a third vortex rotating in the 
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Figura 4 - Streamlines for Stokes’s configuration (Re —» 0), when the cylinder ro¬ 
tates clockwise with (a) c o = 0.1 and (b) uj = 0.5. The rotation of the 
cylinder results in the breakdown of the axial symmetry. As the angular 
speed increases, the lines are deformed, as expected. The wind goes from 
left to right. 


counterclockwise direction begins to form in the region below the horizontal axis, 
see figure 5(c); and this process repeats again and again. 

The method to determine the behavior of the drag force, presented in the 
first section, may now be applied to the case of the rotating cylinder with one 
difference: the volume occupied by the cylinder is initially filled by fluid moving 
exactly as the cylinder, i.e., rotating as a rigid body at time t. Then, starting from 
this configuration, the flow evolves up to t + St. 

In this case, due to rotation, the force on the cylinder has a non zero com¬ 
ponent in the direction transverse to the wind. The components of the force on a 
cylinder that rotates with uj — 0.1 are shown in figure 6, where the forces are expres¬ 
sed in terms of drag and lift coefficients. As the drag coefficient, the lift coefficient 
is an adimensional coefficient widely used in fluid dynamics and defined as (6) 

fi _ Flift 

L - \pV 2 DV 

where Fuf t is the transverse-to-the-wind component of the drag force. The emergence 
of the transverse non-zero component is due to an asymmetry in the boundary layer 
3 separation, triggered first near the bottom of the cylinder (as in figures 4 and 5, 
the cylinder is considered to rotate clockwise). This is a clear manifestation of the 
well known Magnus effect (6). 

As a net result, the von Karman street as a whole is deflected downwards, 
explaining the non null average in figure 6, curve (b), and figure 7, where lift coeffi- 

The boundary-layer results from the adherence of the air molecules to the object, say a cylinder. 
Due to viscosity, the adherence is partially transmitted to the molecules situated farther from 
the cylinder; this defines a region that moves with the cylinder, called boundary-layer or layer 
of Prandtl (6). 


( 6 ) 
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Figura 5 - Starting with Stokes’s laminar configuration for the clockwise rotating 
cylinder, the wind is switched on at t — 0, with a fixed Reynolds number 
Re = 1,000 and oj = 0.5. After a transient time, the velocity field changes 
continually and periodically. Figure (a) shows a vortex that has just 
formed behind the cylinder, rotating in the counterclockwise direction. It 
then moves away to the right. Some time later, another vortex is formed 
behind the cylinder, rotating in the clockwise direction, figure (b). Later, 
figure (c), the velocity field returns to the configuration observed in figure 
(a). The whole process is periodically repeated. At this fine scale close to 
the cylinder back side, the behavior is the same observed in the case of a 
fixed cylinder. However, behind this region there is a long von Karrnan 
street of already previously formed vortices, now bended downwards in 
this slowly rotating case considered here. 


dents as a function of time for various small angular speeds of the rotating cylinder 
are shown. The oscillations are due to the appearance of successive vortices spin¬ 
ning in alternate senses, forming the street. Close to the cylinder, a counterclockwise 
vortex pushes it downwards, subtracting an extra force from the average value, see 
figure 5(a) corresponding to a minimum of the Magnus force in figure 6, curve (b). A 
clockwise vortex, see figure 5(b), pushes the cylinder upwards, adding an extra force 
to the average, corresponding to a maximum of figure 6, curve (b). It is interesting 
to observe that the drag coefficient also oscillates, hardly visible in figure 6, curve 
(a) (see also (5, 25)). 

For larger angular speeds (oj ~ 2.0), it is known that the street of vortices 
tends to disappear (26, 27, 28, 29, 30, 31, 32). Figure 8 shows the behavior of the 
lift, coefficient for some values of uj larger than before. Our results indicate that, 
for angular speeds larger than oo ~ 1.0, the high-frequency oscillations related to 
the sequence of successive vortices are attenuated as the angular speeds increases 
and tend to disappear for angular speed values close to oj ~ 2.0, replaced by much- 
smaller-frequency oscillations. 
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Time 

Figura 6 - The drag force on a rotating cylinder inside the wind tunnel with Rey¬ 
nolds number Re = 1,000, and oj = 0.1. The drag coefficient (a) shows 
the same behavior as before, in the case of fixed cylinder. During the 
transient stage, the lift coefficient (b) also increases, reaches a maximum, 
and decreases. However, this component does not stabilize. Instead, it 
oscillates around a non null mean value. The oscillation is due to the 
alternate formation of counterclockwise and clockwise vortices, being a 
signature of the von Karrnan street. For small values of Re, up to a few 
dozens, both components of the drag force don’t oscillate, a signature of 
the absence of successive vortices. 



Time 


Figura 7 - Lift coefficients as a function of time, for various small angular speeds of 
the rotating cylinder with Re = 1, 000: (a) uj = 0.1; (b) c o = 0.2; (c) uj = 
0.3; (d) u = 0.4; and (e) u = 0.5. After the transient stage, the Magnus 
force oscillates around a non null average value. 
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Time 

Figura 8 - Lift coefficients as functions of time, for various angular speeds of the 
rotating cylinder with Re = 1,000: (a) uo = 1.0; (b) ui = 2.0; and (c) 
co = 2.4. When oo ~ 2.0, the high-frequency oscillations are replaced 
by much-smaller-frequency oscillations, due to the absence of the von 
Karman street. 


4 Final Remarks 

In this paper, we introduce a new method to determine the dynamic behavior 
of the drag force acting on an object inside a wind tunnel. 

Using this method, it was revealed that, in the case of a static cylinder, 
the drag coefficient increases, reaches a maximum value, decreases and stabilizes. 
This behavior agrees with the experimental results presented in various references 
(16, 17, 18, 19, 20, 21, 22, 23). Moreover, the drag forces we found are in excellent 
agreement with known experimental data. These two facts attest to the validity and 
effectiveness of the method. 

In the case of a slowly rotating cylinder, the drag coefficient behaves in the 
same way as in the case of the static cylinder. However, the lift coefficient is not 
null: it grows, reaches a maximum, decreases and fluctuates around a non null value. 
This is a clear manifestation of the well known Magnus effect (6). Moreover, also in 
this case, the method is even capable of identifying a transition scheme where the 
von Karmam street of vortices disappears as the angular speed increases. 

The method, as we hope it has been clear, is simple, because it requires no 
knowledge about the pressure distribution and the shear stresses or the gradient of 
the fluid velocity field on the object surface. No approximation is used (except the 
numerical grid discretization) and no convective correction is necessary. The method 
is also applicable to objects of general shape. 
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6 Appendix 

In this appendix, we demonstrate mathematically the validity of the method 
described in the first section. 

The object (cylinder or any other) inside the wind tunnel is fixed there by 
some mechanical device exerting a force — F on it in order to counterbalance the 
fluid drag force F. At time t one can imagine the object replaced by fluid, also at 
rest. Now, the extra portion of fluid would move. A fluid element d 3 r inside the 
former object volume will suffer a force d /(f) exerted by the rest of fluid around it, 
and thus would acquire a velocity Sv after a small time interval St. Newton’s law 

—* A 7? 

d/(r)=pd 3 r—, (7) 

can be applied, where p is the fluid density. 

I11 order to avoid this extra movement keeping the extra portion of fluid at 
rest, some “magical” device must exert a second force —d f(f) on the fluid element 
d 3 r. Considering the whole extra volume also at rest, one can determine the resultant 
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force this “magical” device must exert on the whole volume V Q b formerly occupied 
by the object 


—F 



( 8 ) 


Now, this “magical” device becomes not magical at all, it is just the same 
already quoted mechanical device fixing the object inside the wind tunnel. Therefore, 
the drag force exerted by the fluid on the object is 


F 


P 


f Vob d 3 r Sv 
St 


(9) 


equation which defines our method. 

Note that the imagined replacement of the object by static fluid is made only 
at the time t, no velocity differences appear. In other words, the above equation is 
exactly the result of Newton’s law as applied to the object. No error, no approxi¬ 
mation. Only one source of inaccuracies may appear: the adopted time interval St 
is not infinitesimal, it must be chosen small enough in the same standards as the 
also finite time interval adopted by the numerical method used to solve the Navier- 
Stokes equation. Therefore, the possible inaccuracies of our method to calculate 
the drag force come exclusively from the numerical method adopted to solve the 
Navier-Stokes equation. These possible inaccuracies are of course controllable. 
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